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ABSTRACT 

In the course of our ongoing search for serendipitous high-redshift Lya emission in 
deep archival Keck spectra, we discovered two very high equivalent width (W^^*^^ > 450 
A, 2a) Lya emission line candidates at z ~ 3 in a moderate dispersion (A/AA ~ 1200) 
spectrogram. Both lines have low velocity dispersions {a^ ~ 60 kms^^) and decon- 
volved radii r ^ 1 h'^^ kpc. We argue that the lines are Lya, and are powered by stellar 
ionization. The surface density of robust, high equivalent width Lya candidates is esti- 
mated to be ~ 3 lb 2 arcmin"^ per unit redshift at z ~ 3, consistent with the estimate of 
Cowie &: Hu (1998). The Lya emission line source characteristics are consistent with the 
galaxies undergoing their first burst of star formation, i.e., with being primeval. Source 
sizes and velocity dispersions are comparable to the theoretical primeval galaxy model 
of Lin & Murray (1992) based on the inside-out, self-similar collapse of an isothermal 
sphere. In this model, star formation among field galaxies is a protracted process. 

Galaxies are thought to be able to display high equivalent widths for only the first ~ 
few X 10'' yr. This time is short in relation to the difference in look back times between 
z = 3 and z = 4, and implies that a substantial fraction of strong line-emitting galaxies 
at z = 3 were formed at redshifts z < 4. We discuss the significance of high-equivalent 
width Lya-emitting galaxies in terms of the emerging picture of the environment, and 
the specific characteristics of primeval galaxy formation at high redshift. 

Subject headings: Cosmology: observations — galaxies: compact — galaxies: evolution 
— galaxies:formation — galaxies: starburst — line: profiles 
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1. Introduction 

Early theoretical models of primeval galaxies suggested that their Lya lines would be highly 
luminous and diffuse. But searches on 4m-class telescopes confirmed no such emission lines (Pritchet 
1994; Thompson &; Djorgovski 1995) and the search for galaxies at large look-back times turned to 
other strategies. One highly successful technique is color-selection of the Lyman break in galaxies 
of redshift z ~ 3 (Stcidcl & Hamilton 1992; Steidel et al. 1996). Lyman-brcak galaxies (LBGs) are 
strongly star- forming, and compact relative to local C* galaxies (Giavalisco et al. 1996b). However, 
LBGs are not generally thought to be primeval, in part because they have lower Lya equivalent 
widths than expected of primeval galaxies, in part because interstellar metallic absorption lines 
(e.g., Si II A1260, Al II A1670) are characteristic of most LBG spectra. 

The theoretical upper limit equivalent width for Lya emission produced by stellar photoion- 
ization is thought to be W^^^ ~ 100 - 200 A (Chariot & Fall 1993), and can be expected to last 
only a few x 10^ yr for a constant star formation rate (SFR). Several examples of high-redshift, high 
equivalent width (non-AGN) Lya emission have been reported (e.g., Dey et al. 1998; Weymann 
et al. 1998; Hu et al. 1999). Few secure examples of isolated, high equivalent widths that are not 
Lya exist in the literature (e.g., Stockton & Ridgway 1998; Stern et al. 2000). About half of LBGs 
have Lya in emission (Steidel et al. 1996, 1998; Lowenthal et al. 1997), with rest equivalent widths 
~ 5 — 35 A. Likewise, about half of local actively star-forming galaxies have Lya in emission, with 
rest equivalent widths in the range Wx 2± 10 — 35 A (Giavalisco et al. 1996a), though a few have 
Wx > 50 A. 

Local studies show that whether Lya is seen in emission or absorption may depend heavily on 
the kinematics of the surrounding neutral H I halo (Legrand et al. 1997; Kunth et al. 1998a) and 
the chance geometry of neutral gas and dust (Giavalisco et al. 1996a). However, VFJ*'^* is shown to 
be weakly-correlated with metallicity (Terlevich et al. 1993; Giavalisco et al. 1996a). The selective 
absorption of resonantly scattered Lya photons, whose path length in the cloud is much longer than 
that of photons of modestly different wavelength, is thought to be the major factor in quenching the 
Lya line. The partially extincted Lya line is characteristically asymmetric (e.g., Dey et al. 1998; 
Kunth et al. 1998a; Legrand et al. 1997), presumably caused by the absorption of photons from 
the blue side of the line by dust in an expanding shell, and compounded by the addition of red- 
shifted backscattering off the neutral shell expanding away from the observer (e.g., Legrand et al. 
1997). Though dust-free galaxies are the most likely to have high equivalent width Lya emission, 
it is clear that a more evolved galaxy can also have high equivalent width Lya emission emerging 
within specific zones of active star formation due to an incomplete covering factor. Such emission 
may appear spatially offset from the region of star formation (e.g., Bunker et al. 2000). Averaged 
over the whole galaxy, the continuum supplied by recent star formation with locally quenched Lya 
emission will reduce the integrated line equivalent width. 

Deep spectra taken at the Keck telescopes regularly reveal serendipitous high-equivalent width, 
isolated emission lines (e.g., Stern et al. 2000). Indeed, the first confirmed galaxy at z > 5 was 
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discovered in this manner (Dey et al. 1998). Steidel et al. (1998) obtained a serendipitous spectrum 

of a strong Lya line with extremely faint continuum near the {z) ~ 3.09 structure in the SSA22 
field. Strong emission lines can also be effectively targeted by using narrow band imaging, in 
conjunction with broad band imaging. Recent observations, using narrow band imaging with a 
"strong equivalent width" {W^^ > 77 A) criterion (Cowie & Hu 1998) and follow-up spectroscopy 
(Hu et al. 1998) have disclosed a population of galaxies with what is thought to be Lya in emission. 
The galaxies tend to have very weak, sometimes undetectable, contimia. They are also are quite 
compact, as noted of the z 2.4 Lya emitters found in HST searches (Pascarelle et al. 1996; 
Pascarelle et al. 1998). Cowie & Hu (1998) estimate a surface density of 3.6 arcmin^^ per unit 
redshift for flux densities j > 2 x 10~^^ ergscm~^s~^, corresponding to line luminosities Ci,ya > 
1.8 X 10^^ ergss"-*- for an Einstein-de Sitter Universe with Hq = 50kms~^ Mpc~^ at z «i 3.4. 

We are conducting a search for serendipitous emission lines in deep archival Keck spectra. 
Such a serendipitous search is an emission-line flux-limited survey, and may find emission lines 
over a wide range of redshifts. Narrow-band surveys are sensitive to candidates over only small 
redshift ranges, and require follow-up spectroscopy to discriminate stellar Lya lines from metal and 

AGN emission lines. The identification of a population of primeval galaxies may provide important 
information about the cpoch(s) of galaxy formation, including data relevant to the integrated global 
star formation rate, and the luminsity function of galaxies to photometric limits fainter than that 
accessible to color-selection surveys. 

One drawback faced by all Lya emission line searches is possible confusion with other isolated, 
high equivalent width lines such as [O II] A3727 and Ha. Experience has shown that [O II] A3727 
is often mis-identified as Lya. When continuum is detected, we may use the Balmer series or the 
continuum depression to discriminate, and when absent we may search for another line, such as 
N V A1240 or [O III] A5007. However, N V A1240 is rarely detected in non-AGN, and the search for 
[O III] A5007 often fails because these lines would either fall out of the range of the spectrograph, or 
in the high-noise infrared part of the spectrum which, combined with the large intrinsic dispersion 
in [O II] A3727/[0 III] A5007 (Kennicutt 1992), can make an unambiguous verification untenable. 
Alternatively, one may use the known tendency of Lya lines to display a P-Cygni profile to identify 
this resonance line. However, asymmetry cannot be a completely necessary criterion for Lya 
emission lines since some local sources apparently have symmetric, or "pure" profiles (Kunth et al. 
1998b). When emission lines are unresolved, and have undetected continua, the last recourse is the 
knowledge that [O II] A3727 rarely has an equivalent width greater than 100 A (but see Stern et al. 
2000). For these reasons, unambiguous Lya lines are rare, and worthy of extra scrutiny. 

We serendipitously found two isolated emission lines on a single slitlet from a mask centered 
on the SSA22 field, andq determined that they are in fact Lya lines — examples of galaxies 
which are more compact and have higher Lya equivalent widths than known LBGs. We report 
on their equivalent widths, intrinsic radii, velocity dispersions, and surf ace/ volume density in or- 
der to place them in the picture of the evolving Universe emerging from studies of high-redshift 
galaxies. The source targeted by the slitlet in which we found our emission line galaxies is the 
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color-selected Lyman-break galaxy SSA22 C17; we analyze it in parallel for comparison. We adopt 
Ho = bOh^o kms~^ Mpc~^, fio = 1; and A = unless otherwise stated. At z = 3, for $7 = 1(0.1), 
I'.'O corresponds to 7.3 (12.3) h'^Q kpc. 

2. Observations and Data Reduction 

We have obtained deep, moderate-dispersion spectra of z ~ 3 LBGs with the aim of making 
detailed studies of the ages, kinematics, dust-content, and abundances of the LBG population 
(Dey et. al., in preparation). The data were taken during the years 1997 to 1999 with the Low 
Resolution Imaging Spectrometer (Oke et al. 1995) at the Cassegrain focus on the Keck II telescope. 
The camera uses a Tek 2048^ CCD detector with a pixel scale (^.'212 pix'^ On UT 1997 September 
10 we obtained moderate-dispersion multislit spectra of LBGs in the SSA22 field, using the 600 
lines mm^-*^ grating blazed at 5000 A. Slitlet widths were 1'.'25, resulting in a spectral resolution 
of ~ 4.4 A (FWHM) for sources filling the slitlet, and a resolution A/AAfwhm ~ 1200. Objects 
not filling the slitlet will have higher spectral resolution. For a spatially unresolved source - one 
blurred only by atmospheric seeing - we estimate a resolution of ~ 2.8 A (FWHM) based on Lya 
absorption systems in the quasar SSA22 D14 also observed on the slitmask reported herein. Four 
of the integrations totaling 7200s were of excellent quality, with seeing along the spatial axis of 
0'.'78. Additional observations were made on UT 1997 September 12 and UT 1999 June 14. Here 
we concentrate on the data of UT 1997 September 10 for which the signal-to-noise ratio (S/N) is 
highest. 

The data reductions were performed using the IRAF package, and followed standard slit spec- 
troscopy procedures. Flat-fielding, sky subtraction, cosmic ray removal, and aperture extractions of 
the slitmask data were facilitated by the home-grown software package BOGUS, created by D. Stern, 
A. J. Bunker Sz S. A. Stanford. Wavelength calibration was performed using a HgNeAr lamp, 
employing telluric sky lines to adjust the wavelength zero-point to the data frames, compensating 
for any drift in the wavelength coverage. Flux calibration was performed using observations of Wolf 
1346 (Massey & Gronwall 1990). 

3. Spectroscopic Results 

We find two serendipitous, isolated emission lines, ser-1 and ser-2, in a 100" long slitlet centered 
on the LBG SSA22 C17 {z = 3.299; see Fig. 1). In lower dispersion spectra, it is often impossible 
to distinguish between high equivalent width forms of [O II] A3727 and Lya, unless there is an 
evident continuum depression, a line asymmetry, or other distinguishing features. Though only 
one of these two serendipitously discovered sources (ser-2) displays evidence of the asymmetry 
characteristic of high-rcdshift Lya emission in LBGs, we argue that Lya is indeed the most likely 
interpretation for both. We simulate the [O II] AA3726, 3729 doublet using our 2.8 A resolution 
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Fig. 1. — Two-dimensional spectrogram of the slitlet in the SSA22 field targeting the LBG 
SSA22 C17 at z = 3.299. Two strong line emitters arc serendipitously identified, ser-1 and ser-2. 
The high-equivalent widths, narrow velocity widths, and lack of secondary spectral features strongly 
argue that these are Lya emitters at z ~ 3. Note the foreground continuum source and the residual 
of the [O I] A 5577 A skyline. 



(for an unresolved source) with 50 kms~^ intrinsic velocity dispersion, introducing the Poisson 
noise appropriate to the observed flux of the median of scr-1 and scr-2 (sec Fig. 2). We conclude 
that [O II] AA3726, 3729 is inconsistent with our lines, as the oxygen doublet would be marginally 
resolved and has a significantly greater FWHM. Further, the absence of associated emission lines 
argues against H/3 and [O III] AA4959, 5007 interpretations. Finally, the lines are shortward of 6563 
A so Ha is not a viable identification. In the following, we assume what is most certainly the case 
— that these are in fact Lya emission lines. Notably, these spectra display no detectable continua. 

The spectra of the serendipitous emission lines ser-1 and scr-2 arc presented in Fig. 3, together 
with that of the LBG SSA22 C17. The redshifts of these galaxies, their emission line wavelengths, 
fluxes and attributed luminosities are given in Table 1. It should be noted that the fluxes and 
luminosities represent only that part of the atmospherically smeared image of the galaxy falling 
within the slitlet, and are in all likelihood under-estimates as these serendipitously-identified sources 



Table 1: Emission Line Properties 



Object 


(A) 


z 


J-17 
(cgs) 


(cgs) 


(A) 


(A) 


C^si) kpc) 


{ km s ^ ) 


ser-1 


4889.3 


3.022 


1.85 


1.29 


> 550 {2a-) 


> 137 (2(t) 


0.8 ±0.6 


62 ± 9 


ser-2 


5296.8 


3.357 


1.35 


1.19 


> 470 {2a) 


> 109 (2(t) 


1.1 ±0.5 


47 ±9 


C17 


5226.4 


3.299 


2.60 


2.20 


35.6 


8.3 


2.8: 


82: 



Notes. — (a) Emission line fluxes are in units of 10 ergs cm ^. (b) The Lya line luminosities axe in units of 
10*^ h^Q ergss^^. (c) radius based on deconvolution of source FWHM. (d) Deconvolved velocity dispersion assumes sources 
fill the 1'.'25 width of slitlets. 
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Fig. 2. — The [O II] AA3726,3729 doublet line profile from a source with a 50 kms~^ velocity 
dispersion, as it would appear in its rest frame, observed with an effective 2.8 A resolution of a 
point-source emitter. Poisson noise has been added commensurate to that observed in the LRIS 
spectra. The observed emission lines of ser-1 and ser-2 have been referred to the same wavelength, 
with their amplitude adjusted to agree with that of the [O II] AA3726, 3729 lines. The line profiles 
of ser-1 and ser-2 are clearly incompatible with an [O II] AA3726, 3729 interpretation. 



arc not necessarily centered in the slitlets. The sources scr-1 and scr-2 arc at rcdshifts z = 3.022, 
and z = 3.357, respectively. We note that ser-2 is at nearly the same redshift as the nearby quasar 
SSA22 D14 (Steidel et al. 1998). 

To date, our serendipitous search has reviewed five slitmasks (~ 0.75 arcmin^) and five longslit 
spectra (~ 0.30 arcmin^). In addition to the two lines presented here, we have found one other 
definite Lya line aX z = 2.946 with a line flux j = 3.32 x 10~^^ ergscm"^ and line luminosity 
jC-Lya = 2.2 X crgss"^. In a total search area slightly over 1 arcmin^, this suggests that the 

surface density of robustly identified high equivalent width Lya emission lines within the redshift 
range, 2.5 < z < 3.5 is very roughly of the order 3 it 2 arcmin^^ per unit redshift. This is consistent 
with the estimate of Cowie & Hu (1998) who find a surface density of « 3.6 Lya-emitters arcmin"^ 
per unit redshift to a higher limiting flux density of 2 x 10~^^ ergscm~^s^^. The full survey 
including candidates at higher redshift will be discussed in a subsequent paper. 
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Fig. 3. — Spectra of two serendipitously identified Lya-emitters in the field of SSA22, with the 
known LBG SSA22 C17. The abscissa on the left refers to the first two plots, while the abscissa on 
the right refers to the right-hand plot. Inserts illustrate the characteristic asymmetric Lya profile 
of C17 and ser-2. 



3.1. Intrinsic Angular Diameters and Velocity Dispersion 

We derive the intrinsic spatial FWHM assuming that the intrinsic, atmospheric, and instru- 
mental FWHM add in quadrature. Based on on the results of multiple measurements, we assume 
the la observational and intrinsic FWHM measurement errors to be ~ 0.25 pix, and run 100 trials 
based on the observed central values. We find deconvolved FWHM for ser-1 and ser-2 of 0'.'24 
lb O'.'lS and O'.'SG it O'.'lS, respectively. At these rcdshifts, a compact galaxy with r = Ih^Q kpc 
would have an intrinsic angular diameter of ~ 0'.'27 / O'.'lS {^l = 1 / 0.1). We note, however, that 
the deconvolved spatial FWHM of C17 at the Lya line, and in the continuum red ward of the line 
are ~ 0^'90 and 0^'40, respectively. The difference is likely due to scattering of Lya photons in the 
H I halo. Hence it is conceivable that the intrinsic sizes of the stellar components of ser-1 and ser-2 
are ~ 50% smaller than the values obtained above, or r < 500h^Q pc. 

We measure the instrumental resolution (for objects filling the slitlet) from the line profile of 

larriD 2 OSO 

lamp spectra. We find FWHMj^^^^ ~ 4.4 A. For spatially unresolved sources, we find FWHMj^^^j, ~ 
2.8 A, measured from the Lya forest of the QSO SSA22 D14. The deconvolved velocity width is 
given by the equation, 

V(FWIIMobJ-^ - (FWIIM„,,trF 
2.354 Aobs ' 

where c is the velocity of light. We find that is consistent with ~ 55 it 8.5 kms~^ for both 
serendipitous objects assuming the sources fill the slitlet (see Table 1). For unresolved sources, 
perhaps more appropriate for the small angular extent of these sources, ^ 95 kms^^. Bear in 
mind, however, that these values refer to the kinematics of the circum-galactic H I cloud, rather 
than stellar velocity dispersion, which will be less (see related discussion in §5). 
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3.2. Lya Equivalent Widths 

The measured continua for the serendipitously-identified Lya-emitters were consistent with 
zero counts, confirming the visual impression of ser-1 and ser-2 seen in Figs. 1 and 3. Therefore, 
direct spectroscopic measure of the hne equivalent widths is not possible. We can, however, deter- 
mine lower limits on the equivalent widths using the la Tinccrtainty of the continua. The Poisson 
noise is the square root of the photon counts per pixel, and wc denote its average per pixel by 
cTpix- The summed noise in a aperture of width w and length I is api^Vwl. We have I different 
approximations of the value of the continuum, so the noise in the continuum cTcont is, 

•^cont ~ f^pix 

However, w is not a free parameter, but must be chosen so as to maximize the S/N. Since our 
sources are small, we assume their undetected continua have a Gaussian shape along a spatial cut, 
and that their FWHM is the same as the seeing (~ CK'TS). For the Lycc line, this would imply a 
source of zero intrinsic extension. But since the FWHM of the continuum is probably significantly 
less than that of the Lya line, as we found with SSA22 C17 (§3.1), it may in fact be a slight 
over-estimation of the continuum FWHM. We find the S/N ratio is maximized when w = 1.165x 
FWHM. For / we chose 100 pixels 126 A). Then, 

^obs^2a)>-^. (3) 

cont 

For the 1 a continuum uncertainties, CTcont — 1.7(1.4) x 10~^° ergscm"^ s~^ A ^ for ser-1 (ser-2). 

We find the 2a lower-limits on the observed equivalent widths are W^^^ > 550 A (ser-1), and > 470 
A (ser-2). The rest-frame values are tabulated in Table 1. Considering the apparent absence of N V 
A1240 and C IV A1549 and the low ISM velocity dispersion of these sources, the data is consistent 
with stellar sources, and only marginally with a possible weak AGN. We discuss this further in §4. 
SSA22 C17 has measured VF^^^ ~ 35.6 A. 

3.3. Star Formation Rates 

Assuming the Lya line flux is powered by star formation, the strength of this line can be used 
to infer a star formation rate for a young galaxy. We first estimate the Ha luminosity corresponding 

to the measured Lya luminsity: the low-density case B ratio Lya/H/? is approximately 25 in the 
absence of dust (Ferland & Osterbrock 1985), while Ha/H/3 ~ 2.8 for hot stars (T ~ 1.2 x 10"^ K; 
Osterbrock 1989). Kennicutt (1983) provides the conversion from Ha luminosity to the SFR, 
resulting in a SFR - CLya relation of M {Mq yr^^) ~ lO'^^^^Lya ( 

ergss ). This results in SFRs 

of ~ 1.3/1.2 Mq yr~^ for ser-l/ser-2, respectively. Because of the possibility of dust, or incomplete 
absorption of ionizing radiation by the gaseous halo, these numbers are lower limits. We believe slit 
losses are minimal as the lines are nearly unresolved and do not fill the slit. Line fluxes from two 
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other nights of observations support this conclusion. The continuum flux of the LBG SSA22 C17 
at 1500 A was used to determine the SFR of M = 27.7 Mq yr'^ (Madau et al. 1998), based on 
a Salpeter IMF with mass range 0.1 < M/Mq < 125. However, based on the synthetic spectra of 
Leithercr ct al. (1995), with a continuous SFR at ~ 10 Myr of age, the SFR is ~ 1.6 times larger, or 
~ 45 Mq yr~^. We expect that the true value is somewhere in between. The Lya line luminosity 
of C17 (Table 1) would imply a SFR of 2.2 M© yr. 

3.4. Projected Photometric Magnitudes 

We use the spectral slope of a 10^ year old starburst with a rest equivalent width of 150 A(Charlot 
& Fall 1993) to calculate the expected continua of ser-1 and ser-2. We would thus expect continuum 
magnitudes oi V = 27.8/28.0 for ser-l/ser-2, ox V = 27 A/27. 7 when the Lya line is in included. 
The 2(7 upper limits on the continua of ser-1 and ser-2 (§3.2) are only ~ 30% greater than the 
continua implied by W^f^* = 150 A. 

4. Discussion 

What is the physical origin of these isolated, high equivalent width emission lines? N V A1240 
or C IV A1549 are not found in the spectra of ser-1 or ser-2. The ratio of the high-ionization line 
flux relative to Lya is < 1/15 and < 1/20 for N V A1240 and C IV A1549, respectively, significantly 
less than ratios measured from composite AGN spectra (e.g., Stern ct al. 1999). The small velocity 
dispersions of ser-1 and ser-2 ((j„ ~ eOkms^i) and the lack of associated N V A1240 or C IV A1549 
emission argue against the presence of an AGN, while the high surface brightness is inconsistent 
with photoionization by the metagalactic flux (e.g.. Bunker et al. 1998). Thus we assume the lines 
owe their existence to stellar sources. Below we argue that the galaxies are primeval, and attempt 
to fit them into a cogent picture of the evolving Universe at 2; 2± 3. 

4.1. The Production of High Equivalent Width Lya Emission 

The production of a strong Lya line requires two things: first there must be a large column 
of neutral hydrogen which can intercept almost all photons with A < 912 A, and second, Lya 
radiation must escape the cloud. The large cross section of H I for Lya photons makes Lya 
extremely vulnerable to dust quenching. We consider the scenario of Tcnorio-Tagle ct al. (1999) 
in some detail as one model of the early phases of a starbursting protogalaxy. Tenorio-Tagle 
et al. (1999) suggest that after an initial stellar ionization of dynamically quasi-relaxed gas and 
attendant high equivalent width "pure" Lya emission (i.e., a symmetric line profile), enriched stellar 
winds sweep up gas. The compressed gas recombines, forming a shell of neutral gas which absorbs 
emerging Lyman continuum radiation. Lya radiation is then emitted and scattered within the shell. 
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The expansion of the compressed shell of gas results in the preferential scattering, and probable 
absorption, of the blue-wing of the Lya line. The LBG SSA22 C17 (Figs. 1, 3) is a good example 
of this common feature of LBGs; in the insert we see a trough with center ~ 350 kms~^ blueward 
of the line center, while at ~ 650 kms^^ we see a second maximum. Here, perhaps, is palpable 
proof of a strong wind and a neutral shell, with the blueward emission line probably caused by the 
recombination of shocked gas associated with the expanding shell (Tenorio-Tagle et al. 1999). In 
this scenario, the time scale for high equivalent width Lya line emission must be only modestly 
greater than the lifetime of the most massive stars (at most ~ 10^ yr). 

This model is based upon a single, central starburst driving the Lya emission — an assumption 
which may not generally hold in real galaxies. High-resolution VLA/optical observations of five local 
star- forming blue compact dwarfs reveal a history of spatially distinct star forming regions in various 
stages of evolution (van Zee et al. 1998). Perhaps new star- forming regions in the same galaxy could 
display high equivalent width Lya emission. This effect appears to be seen in a z = 4.04 lensed 
galaxy whose star formation regions are subjected to a detailed lensing chromatography study by 
Bunker et al. (2000) . Their color analysis shows a sequence of non-coeval star formation sites with 
ages ranging from zero to ~ 100 Myr. "Nodes" of star formation, showing Lya in absorption, 
have spectral energy distributions consistent with a burst age of 20 — 30 Myr. There is, however, 
an outlying region of very high equivalent width Lya emission (Wf > 100 A, 3a) thought to 
be resonantly-scattered Lya photons escaping from the adjacent star- forming node. Yet summed 
over the whole galaxy (as we would see it without the benefit of gravitational lensing), the Lya 
equivalent width is modest. 

We therefore suggest that, for apertures containing entire young and compact galaxies, high 
equivalent width stellar Lya emission {WJ^^^ > 100 A) indicates a primeval or near-primeval galaxy. 
We conclude that it is very likely that ser-1 and ser-2 are primeval galaxies — by which we mean 
they are undergoing their first significant burst of star formation. 

4.2. The Surface Density of High-Redshift Lya-Emitters 

How do ser-1 and ser-2 fit within the context of surveys of star forming galaxies at high 
redshift? The luminosity function (LF) of galaxies at high redshift is a useful comparative tool for 
our purposes. Steidel et al. (1999) find a steep faint-end slope (a ~ —1.6) for the observed 7?.-band 
LF of star-forming galaxies at (z) ~ 3.04, implying that a large fraction of the UV luminosity density 
is produced by galaxies fainter than the spectroscopic limits of current Earth-based photometric 
surveys. By virtue of the association of UV luminosity with SFR, this implies that a large fraction 
of the star formation history is unobserved. When integrated to the apparent magnitudes expected 
of the continua our serendipitous emission- line galaxies (JZ ~ 28 mag) , their UV luminosity function 
implies number densities ~ 20 times, and integrated star formation density ~ 4 times, higher than 
the classical LBGs integrated to the completeness limits of their ground-based survey, TZ = 25 mag. 
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As we have noted, because the Lya emission hnes of such sources as our serendipitous objects 
are Ukely to be observable only for a very small fraction of 1 Gyr, the modest numbers of high 
equivalent width Lyct sources are tangible evidence for this many-fold larger population of slightly 
older, faint continuum, star forming galaxies with self-absorbed Lyo inaccessible to even the deepest 
spectroscopic or color selection surveys. If we are to believe that high equivalent width Lya can 
only be exhibited for the first ~ 5 x 10^ yr of a galaxy's life, and that their surface density is 
~ 3 arcmin"^ per unit redshift, then if their distribution is uniform between 3 < z < 4, there must 
be ~ 30 (50) galaxies (fi = 1 (0.1)) per square arcminute formed in that interval. This follows 
because there are 10 (17) intervals of 5 x 10^ yr between z = 3 and z = A. The surface density of 
newly formed galaxies would thus be ~ 30 times larger than the LBG surface density. One must 
integrate the UVLF to TZ = 28.5 to achieve that surface density. Thus, available surface density 
estimates of primeval (high equivalent width Lya emission) galaxies are consistent with the claim 
that most galaxies at z ~ 3 were formed at redshifts less than z = 4. However, the integration 
of the z ~ 4 LF (Steidel et al. 1999) over luminosity limits identical to that at z = 3 shows that 
the implied co-moving number densities increase only slightly for 0, in the range of 0.1 < < 1.0, 
apparently consistent with very low net galaxy formation. These two lines of research could be 
consistent if the merger rate was nearly equal to the galaxy formation rate. There is evidence that 
the epoch-dependent merger rate was higher in the past by (l + z)"^ with m ~ 1.5 — 3.0 (Windhorst 
et al. 1999), and that mergers were more frequent for the sub-£* galaxies. However, determining 
the size-dependent merger rate that would be needed to leave the shape of the luminosity function 
unchanged while producing only a very modest increase in comoving galaxy density, is beyond the 
scope of this paper. 



4.3. Environment 

The environment of high equivalent width Lya-emitters is of considerable interest. In their 
study of Lyman-break candidates in the SSA22 field, Steidel et al. (1998) display convincing ev- 
idence for redshift "spikes" at {z) = 3.09, and (z) = 3.35. Ser-1 lies ~ 5000 km s~^ from the 
former enhancement, a structure which they suggest is a proto-Abell cluster. The latter redshift 
spike, which is substantially less massive than the former, contains a quasar (SSA22 D14). Ser-2 
is only ~ 190(320)/i5Q^ kpc (projected), for = 1 (0.1), and < 200kms~^ distant from the quasar, 
while both are on order 400 kms~^ beyond the apparent center of the redshift spike^. We note 
that the serendipitously-identified high equivalent width Lya line, SSA22-S1 (Steidel et al. 1998) 
at z = 3.100, is ~ 750kms~^ from the (z) = 3.09 enhancement. Thus, this limited amount of 
data suggests that high equivalent width Lya emitters tend to be loosely associated with density 
enhancements some of which may be protoclusters. Similar clustering of Lya emitters is seen in 



''We modeled the effect of the quasar on ser-2 by extending the observed continuum blueward of the Lya line. 
The extrapolated photon flux is such that its excitation of gas around ser-2 could account for at most ~ 3% of the 
Lya emission of ser-2. 
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studies at ^ ~ 2.4 (Keel et al. 1999). 



5. A Primeval Galaxy Model 

In a speculative vein, we now consider what we have learned of primeval galaxies in the light 
of observations, and the modeling of associated phenomena. Our purpose is to emphasize the 
significance of the above results to the wider field of cosmology. Our study suggests the existence 
of a high volume density of tiny galaxies with high specific SFRs, which may cumulatively have 
rivaled that of the LBGs. We seek here to provide a galaxy formation scenario in which this is 
plausible, and to place this scenario in the larger context of the evolution of the universe. 

In response to the very small effective radii of faint galaxies imaged by HST (e.g., Giavalisco 
et al. 1996b; Pascarelle et al. 1998), it has been shown that an "inside-out" galaxy formation 
scenario may produce galaxies with small scale-lengths (Bouwens et al. 1997; Cayon et al. 1996). 

A theoretical framework for this hypothesis is presented in a paper by Lin & Murray (1992), which 
seeks to model primeval galaxies using two alternative hypotheses of protogalactic clouds. In these 
models, robust self-regulating star formation occurs at a rate inversely proportional to the cooling 
time for hydrogen number densities greater than Uc ~ 4cm~^. Their model A is a homogeneous 

1/3 

cloud which collapses uniformly; the critical density is first reached at a size Rc ^ 6M^{ kpc, where 
Mil is the galaxy mass in units of 10^^ Mq, creating a star-forming region whose size is inconsistent 

with observations of either ser-1 or ser-2. Model B is an isothermal cloud with a core of 500 pc, 
which collapses in an inside-out manner according to the Shu (1977) self-similar collapse model, and 

1/2 —1/2 

reaches the density ric at Rc ^ lAfj^i -^^loo ' where i2ioo is the outer radius of the cloud in units 
of 100 kpc. Model B predicts sizes that are very close to the constrained sizes of our serendipitous 
emission line regions, and makes their comparison with the model appear quite promising. It 
also predicts that protogalaxics will have a relatively constant bolometric luminosity for a period 
of ~ 1.7 Gyr. In their published simulation, the bolometric luminosity log£ ~ 43.5 (erg s^^), 
though the scale-free nature of isothermal spheres could be expected to allow the somewhat more 
modest luminosities observed here. We shall see that this luminosity is consistent with a SFR of 
~ 5.7 Moyr-^ 

An isothermal cloud which is bounded by an external medium of pressure Pgxt ™ay achieve a hy- 
drostatic equilibrium if the ratio of the external to the internal pressure is less than 14.3 (Shu 1977). 
This cloud, known as a Bonnor-Ebert sphere, has a just-critical mass Merit - 0.8a4G-3/2p-y2^ 
where a is the sound speed. This mass may be consistent with galactic masses when the intergalac- 
tic pressure is modest (i.e., Pcxt ^ 1 x 10~^^dyne cm~^). Such a sphere would become super-critical, 
and collapse from the inside out if it were part of a general over-density which had pulled away 
from the Hubble fiow and begun to collapse. Linear theory shows that the turn-around radius 
moves outward at a rate of ~half the Hubble flow (Davis et al. 1980), so that we could expect 
the location of high equivalent width Lya emission to be somewhat within the turn-around radius, 
and to move outward with time. A more extensive discussion of this subject will be reserved for a 
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future publication. 

In order to link the model with observations, we made a rough simulation of a galaxy under- 
going a constant SFR over a period of up to 1 Gyr. We used a Salpeter IMF (cc = —2.35) with 
lower mass limit of 0.1 Mq, and main sequence lifetimes and luminosities (Tables 3-6 and 3-9) 
from Mihalas & Binncy (1981), modeling them as black bodies, which is sufficiently accurate for 
metal- free, non-evolved stars. Unless otherwise stated, our upper mass limit is 125 Mq, and we 
assume that the SFR is 1.2 Mq yr~^. The Lya photon emission rate of a primeval galaxy is taken 
to be 2/3 of the stellar emission of photons more energetic than 13.6 eV, according to the case B 
recombination (Spitzer 1978). We assume that 80% of all newly introduced gas (i.e., over and above 
the ~ 10^ Mq of gas of density ric within the central 500 pc) is transformed into stars immediately. 

The continuum measured in the 75 A region around the Lya emission line, together with 

the Lya luminosity, is used to calculate the equivalent width. We tested the derived equivalent 
widths against the dust-free calculations of Chariot & Fall (1993) for continuous star formation 
and upper mass limits of 80 and 120 Mq, and found values that were quite consistent with values 
attained by interpolating between the bracketing Salpeter slopes which they employed. We find 
that for a 10 Myr old galaxy, = 160/200/250A respectively for 80/100/120 Mq upper limits. 

These derived values drop to ~ 125/155/200A by 100 Myr. The relation of the resulting Lya 
line luminosity to the SFR is in agreement with the case B predictions. The model predicts a 
continuum in the observed visual band of ~ 1.3 x 10"^'' ergscm"^ s~^ A for a z = 3 galaxy with 
M = 1.2 Mq yr-^ This is within the la non-detection limits of the continua of ser-1 and ser-2 
(1.7/1.4 X 10-20 ergscm-^s-i A~^). 

For an upper mass limit of 100 Mq, and M = 1.0 Mq yr~^, the bolometric luminosity at 
10 Myr is ~ 5.5 x 10^^ ergss"^. Thus, the luminosity of the Lin & Murray (1992) model B 
corresponds to M 2± 5.7 Mq yr"^. We calculate the mass to bolometric luminosity ratio of our 
model by adding the mass in stars and dissipated gas, dividing it by the luminosity. The dark 
matter does not participate in the cloud's collapse since it is not pressure supported. It will be 
compressed somewhat by the gradually increasing concentration of baryons in the center, however, 
on the whole, the DM halo should have little effect on the velocity dispersion of stars formed out 
of dissipated gas in the first few tens of Myr since they are formed well within the core of the 
isothermal sphere. We find that the calculated A^/>C 2± 0.01, 0.05, and 0.37 for logt = 7, 8, and 9 
(Gyr) respectively. We note that compact narrow emission line galaxies, which appear to be likely 
lower-redshift counterparts, have mass-to-luminosity ratios, M/jC 0.1 (Phillips et al. 1997). 

The expected stellar velocity dispersion, ag = GMstars+gas/f, implies dispersions of model 
galaxies of from ~ 14 / 37 / , and 110 kms"^ for logt = 7, 8, and 9 (yr), assuming a radius of 500 
pc. The low o"s for the 10 Myr old galaxy is actually reasonable, in view of the finding that the very 
large cluster in the blue compact dwarf NGC 1705 (Ho &: Filippenko 1996) has cr^ ~ llkms"^. A 
low as is also not in conflict with our deconvolved velocity dispersions of ay ~ 60 kms"^ (§3-1) 
since the latter is a measure of both the turbulence of gas (from infall and stellar winds) which give 
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rise to the Lya resonance line, as well as the possibly spatially resolved structure within the slitlet. 
On the basis of the observed and deconvolved radius r of the LBG C17 (Table 1), we find that 
the galaxy mass is Mlbg = 4.3 x 10^ Mq. This galaxy would take ~1.2 Gyr to form at the SFR 
derived for ser-1 or ser-2. However, between z = 3.5 and z = 3.0, there are only 0.265/0.473 Gyr 
for J7 = 1.0/0.1. Thus, LBGs would not be formed in a reasonable amount of time from low-SFR 
objects if SFRs remain constant. Below, we consider the possibility that the SFR may increase 
with time 

Our simple galaxy model has assumed that M is constant. However, a realistic isothermal 
sphere will have a baryonic core, assumed here to be ~ 500 pc. In this case case the infall rate 
will not be constant in time, but given a sound speed of 10 kms~^, may increase during the first 
~ 5 X 10'' yr. In addition, the negligible metallicity expected of primordial gas will make cooling 
times, and the Jeans mass, larger^ during the first few tens of Myr of a galaxy's life. Thus the effect 
of metals ejected from evolved stars in aiding cooling, and of dust in catalyzing the production of 
molecules, is to compound the effect of increasing infall; the SFR and its efficiency can be expected 
to rise substantially by a few tens of Myr. 

On the other hand, recall that the high equivalent width of the Lya emission line of primeval 
galaxies is expected to be quenched over essentially the same time scale. This means that the 
nature, and fate of these pre-LBG/post-primeval galaxies will be very difficult to learn. With an 
ernhanced SFR, the post-primeval galaxies could plausibly grow by a factor of 10-40 into LBGs 
in a few tenths of a Gyr or, by virtue of competition from neighbors, have their growth truncated 
while still in the dwarf stage, or even be cannibalized; environment may have a dominant role to 
play. 

While the mass and SFRs of our purportedly primeval galaxies are much smaller than that 
of LBGs, the specific SFR of the former is on average 10 or 20 times greater than the latter. 
The large numbers of small, young galaxies which are thought to have been formed (§4.2) in the 

interval 3 ^ z ^ 4 might thus be expected to provide a substantial fraction of the star formation 
rate density. In fact, the steep faint end slope of the z ~ 3 UV luminosity function of Steidel et al. 
(1999) implies that up to 40% of the SFR density may be supplied by galaxies fainter than TZ = 27 
{C ^ 0.1) but brighter than the projected continuum magnitudes of our serendipitous sources. 

In addition to the formation of a galactic bulge according to the Lin-Murray process, it is 
expected to be accompanied by the formation of the halo, as infalling clumps of gas dissipatively 
interact with the already accreted gas (Binney 1976; Manning 1999). 



^It has been suggested (Lin & Murray 1992) that the IMF of primeval galaxies may be heavily weighted in 
high-mass stars, which could result in rest Lya equivalent widths significantly larger than 150 A. The truth of this 
contention could be tested if spectroscopic observations deep enough to detect the continua were made. 
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6. Conclusions 

We draw the following conclusions from this study: 

• We suggest that the combination of high equivalent widths, low velocity dispersions, and 
small intrinsic sizes seen in the Lya lines of these two galaxies can be well-explained as 
their being extremely low-metallicity galaxies undergoing an initial burst of star formation — 
hence, primeval galaxies. These characteristics, markedly different from earlier expectations, 
suggests an inside-out, rather than a monolithic collapse formation mechanism. 

• Our data, which imply that most galaxies at z = 3 may have been formed at z < 4, stand 

in apparent conflict with the z = 3 and 4 UVLFs (Steidel et al. 1999), which suggest that 
the comoving galactic density remained relatively constant during this interval. These may 
be reconciled, however, with a large merger rate that preserves the faint-end slope of the 
luminosity function. 

• The agreement of the angular size, and the projected total luminosity of these emission line 
galaxies, with the predictions of the inside-out collapse of a nearly-isothermal H 1 halo (Lin 
Sz Murray 1992) suggests that the latter is a promising basis for modeling galaxy formation. 
Our discussion has suggested that a rising star formation rate with a time-scale of ~ 5 x 10^ 
yr is plausible for many post-primeval galaxies. 

• The isothermal collapse model requires that protogalactic clouds form in relatively isolated 
regions. However, their collapse may be stimulated by the gradually increasing pressures of 
a protoclustcr environment, thus accounting for the apparent weak clustering on large scales 
noted of high-, as well as low-redshift Lya emission line galaxies. 

We thank C. Leitherer, and C. C. Steidel, for useful discussions. We are grateful to our 
referee. Dr. R. A. Windhorst, for valuable suggestions. We also thank A. Dey for assistance with 
the observations. We are grateful for the support of NSF grant AST 95-28536. 

References 

Binney, J. 1976, MNRAS, 177, 19 

Bouwens, R. J., Cayon, L., & Silk, J. 1997, ApJ, 489, L21 
Bunker, A. J., Marleau, F. R., & Graham, J. R. 1998, AJ, 116, 2086 
Bunker, A. J., Moustakas, L. A., &: Davis, M. 2000, ApJ, in press 
Cayon, L., Silk, J., & Chariot, S. 1996, ApJ, 467, L53 



-16- 



Charlot, S. & Fall, S. M. 1993, ApJ, 451, 580 
Cowie, L. & Hu, E. M. 1998, AJ, 115, 1319 

Davis, M., Tonry, J., Huchra, J., & Latham, D. W. 1980, ApJ, 238, L113 

Dey, A., Spinrad, H., Stern, D., Graham, J. R., & Chaffee, F. 1998, ApJ, 498, L93 

Ferland, G. J. & Osterbrock, D. E. 1985, ApJ, 289, 105 

Giavalisco, M., Koratkar, A., & Calzetti, D. 1996a, ApJ, 466, 831 

Giavalisco, M., Steidel, C. C., & Macchetto, D. 1996b, ApJ, 470, 189 

Ho, L. C. & Filippenko, A. V. 1996, ApJ, 472, 600 

Hu, E. M., Cowie, L. L., & McMahon, R. G. 1998, ApJ, 502, 99 

Hu, E. M., McMahon, R. G., & Cowie, L. L. 1999, ApJ, 522, 9 

Keel, W. C, Cohen, S. H., Windhorst, R. A., & Waddington, 1. 1999, AJ, 118, 2547 

Kennicutt, R. 1983, ApJ, 282, 54 

— . 1992, ApJ, 388, 310 

Kunth, D., Mas-Hesse, J. M., Terlevich, E., Terlevich, R., & Fall, S. M. 1998a, A&A, 334, 11 

Kunth, D., Terlevich, E., Terlevich, R., k, Tenorio-Tagle, G. 1998b, in Dwarf Galaxies and Cos- 
mology, ed. T. X. Thuan (Gif-sur-Yvette: Editions Frontieres), in press 

Legrand, F., Kunth, D. F., Mas-Hcssc, J. M., k Lcqucus, J. 1997, A&A, 326, 929 

Leitherer, C, Carmelle, R., & Heckman, T. M. 1995, ApJS, 99, 173 

Lin, D. N. C. & Murray, S. D. 1992, ApJ, 394, 523 

Lowenthal, J. D. et al. 1997, ApJ, 481, 673 

Madau, R, Pozzetti, L., & Dickinson, M. E. 1998, ApJ, 498, 106 

Manning, C. 1999, ApJ, 518, 226 

Massey, P. & Gronwall, C. 1990, ApJ, 358, 344 

Mihalas, D. k. Binney, J. 1981, Galactic Astronomy (San Francisco: W. H. Freeman and Company) 
Oke, J. B. et al. 1995, PASP, 107, 375 

Osterbrock, D. 1989, The Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Mill Valley: 
Univ. Science Books) 



-17- 



Pascarelle, S. M., Windhorst, R. A., & Keel, W. C. 1998, AJ, 116, 2659 

Pascarelle, S. M., Windhorst, R. A., Keel, W. C, & Odewahn, S. C. 1996, Nature, 383, 45 

Phillips, A. C, Guzman, R., Gallego, J., Koo, D. C, Lowenthal, J. D., Vogt, N. P., Faber, S. M., 
& Illingworth, G. D. 1997, ApJ, 489, 543 

Pritchet, C. J. 1994, PASP, 106, 1052 

Shu, F. H. 1977, ApJ, 214, 488 

Spitzer, L. 1978, Physical Processes in the Interstellar Medium (New Work: Wiley) 

Steidel, C. S., Adelberger, K. L., Dickinson, M., Giavalisco, M., Pettini, M., & Kellogg, M. 1998, 
ApJ, 492, 428 

Steidel, C. S., Adelberger, K. L., Giavalisco, M., Dickinson, M., & Pettini, M. 1999, ApJ, 519, 1 
Steidel, C. S., Giavalisco, M., Dickinson, M., & Adelberger, K. L. 1996, AJ, 112, 352 
Steidel, C. S. & Hamilton, D. 1992, AJ, 104, 941 
Stern, D., Bunker, A. J., Spinrad, H., Sz Dey, A. 2000, ApJ, in press 

Stern, D., Dey, A., Spinrad, H., Maxfield, L. M., Dickinson, M. E., Schlegel, D., & Gonzalez, R. A. 
1999, AJ, 117, 1122 

Stockton, A. & Ridgway, S. E. 1998, AJ, 115, 1340 

Tenorio-Tagle, G., Silich, S. A., Kunth, D., Terlevich, E., & Terlevich, R. 1999, MNRAS, 309, 332 

Terlevich, E., Diaz, A. L, Terlevich, R., & Garcia, M. L. 1993, MNRAS, 260, 3 

Thompson, D. & Djorgovski, S. G. 1995, AJ, 110, 982 

van Zee, L., Skillman, E. D., k Salzer, J. J. 1998, AJ, 110, 963 

Weymann, R., Stern, D., Bunker, A. J., Spinrad, H., Chaffee, F., Thompson, R., & Storrie- 
Lombardi, L. 1998, ApJ, 505, L95 

Windhorst, R. A., Cohen, S. H., &; Waddington, I. 1999, in After the Dark Ages: When Galaxies 
were Young, ed. S. S. Holt & E. P. Smith, Vol. 470 (New York: AIP), 202 



